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Abstract 

We study a gauged Chern-Simons-CP(l) system. We show that contrary to previous 
claims the model in the absences of a potential term cannot support finite size soliton 
solution in R 2 . 
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The CP(n) sigma model have been investigated in detail since the early 70's mainly 
as toy models to explore the strong coupling effects of QCD and as effective models of 
some condensed matter systems. In the 2 + 1 dimensional case they have been argued 
to probably play a role in the description of High T c superconductors pp. 

An important issue related to this type of models concern the existence of soliton 
type solutions. For the simplest CP{1) mode topological solutions have been shown 
to exist[2]. Nevertheless, the solutions are of arbitrary size due to scale invariance. As 
argued originally by Dzyaloshinsky, Polyakov and Wiegmann[3j a Chern-Simons term 
can naturally arise in this type of models and the presence of a dimensional parameter 
could play some role stabilizing the soliton solutions. A first detailed consideration 
of this problem was done in Ref.[l] where a perturbative analysis around the scale 
invariant solutions (i.e no Chern Simons coupling k = 0) showed that the solutions 
were pushed to infinite size. Still, it remained open the possibility of non perturbative 
solutions. This issue was considered in Ref.[5] where by direct numerical analysis the 
authors claimed the existence of stable solitons. 

The aim of this note is to reconsider this problem in more detail and as a result 
of our investigation we will show that contrary to the claim made in Ref.[5] that this 
solutions do not exist in R 2 . Still, we will show that the problem can be considered in 
a finite two dimensional disk, which could be a possible scenario in condensed matter 
applications (see [6] for the analogous problem of vortices in Ginzburg Landau theories) 

We begin by considering a (2 + 1) dimensional Chern-Simons model coupled to a 
complex two component field n(x) described by the action 



The subindex D indicates that the region of integration is a disc D of radius R. Here 
D ^ = + iA^ (fi = 0, 1,2) is the covariant derivative and S cs is the Chern-Simons 
action given by 



S 




(1) 




(2) 
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where 

F» v = A u - d„A„ (3) 

The signature of the metric g^ u is (1, —1, —1) and the two component field n(x) is 
subject to the constraint n^n = 1 

Variation of the action yields the field equations 



DpD^n = {n ] DpD»n)n (4) 

^v P F vp = -Jp = i[rJD„n - n{D^\ (5) 
The time component of last equation 

2kF 12 = -J (6) 

is Gauss's law of Chern-Simons dynamics. Integration over the disc gives a relation 
between charge Q = J D d 2 xJo and magnetic flux = f D F\2(i 2 x [7]: 

• — ifl (7) 

Defining the stress tensor as T^ u = and using the equations of motion, the energy 
functional for a static field configuration can be expressed as 

E = J d 2 x(n 2 B 2 + |An| 2 ) , i = l,2 (8) 

Let us consider the following ansatz for the iV soliton solutions [5]: 



n(M=\ y \ |, Mr) = <r) , A r = (9) 

sin(-Tp) 



2 

Using this ansatz Eq.(|8]) becomes 



E = 2 ^j o rdr{n 2 + d r a{r)^ + ^(d r 6(r)) 2 

( N 2 2Na(r)\ 9 Mr), 9 , A 



3 



We shall impose the following boundary conditions at the origin 

lim#(r) = 7r, lim a(r) = (11) 

These conditions imply regularity of the fields at the origin. On the other hand the 
conditions at the boundary of the disk are in principle more general. If we were working 
on the infinite plane, then the natural boundary conditions ensuring finite energy would 
be, 



N 

lim <9(r) = , lim a(r) = (12) 

r— s-oo r->oo r 

For a finite disk it seems reasonable to impose 

N 

lim a (r) = -— (13) 

r^>R K 

which implies quantization of the magnetic flux 

$ = 2tt [ R rdr ^ ^ r)) = -2irN (14) 
Jo r 

For the 6 field we could impose a Dirichlet type boundary condition, 

lim 0(r) = 9(R) (15) 

or a Neumann type boundary condition 

lim ^ = *'(*) (16) 
r^R dr 

The condition 0(R) = is special in the sense that it implies the vanishing of the 
current at the boundary, 



Jrh 



Using the equations of motion, it can be shown that, 



2k 2 h 

where the CP(1) charge density is defined as 



VJS-^Jb + ^jf"" (18) 



jCPil) = » e «( An )t(i>J n ) (19) 

Then, 

B = 2ttJq P{1) - kV.E (20) 

Using that, 

2kV.E = -(3iJ 2 - 8b JO (21) 

$=/ Bd 2 i = 2i / J CP W d 2 x + - / Jdx i (22) 
in io 2 J dD 

We clearly see that the relation between magnetic flux and CP(1) charge involves a 

boundary term related to the current. 

Using the Bogomol'nyi identity [H] |A™| 2 = \{D\ ± iD2)n\ 2 T B ± ^-d%Jj , the 

energy ([8]) becomes 

E = f d 2 x (n 2 B 2 + |(Z>i ± iD 2 )n\ 2 ) T 2vrQ CP(1) (23) 

JDR V J 



I DR 

We see that the energy is bounded bellow by the CP(1) charge. 



E > 2tt / d 2 x Jcp(i) = 2vr|Q 

JDR 



CP(1)\ 



(24) 



This bound is saturated by fields satisfying the first-order Bogomol'nyi self-duality 
equations [8]. 

\(Dx ±iD 2 )n\ 2 = 0, B = (25) 
The field equations corresponding to the ansatz Eq.Q read 



d 2 T a{r) + ^ " ^ " ^ = - 2 (^)4 (26) 



rd r (rd r 6>(r)) + TV 2 sin(0(r)) = -2Nra{r) sin(0(r)) 
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(27) 



If the solutions of ([26]) and ([27]) exist their scale must be set by the quantity k. Following 
Ref.[5], we introduce the dimensionless quantities 



A = na , s = ~ (28) 

K 



in terms of which (1261) and (1271) become 



v> , d s A A . 2/ 9,N 

d 2 A + ^ A = cos 2 (-- 29 

s s z 2 s 



sd s (sd s 9) + N 2 sm(9) = -2NsAsm(9) (30) 
The energy functional (|10p in terms of these new variables reads as 



2 

s^ 2 



E(S) = 2irJ^ sds((~+d s A\ + ~(d t 
( N 2 2NA\ 2 ,9, , 2 \ 

For the origin we choose the following boundary conditions, 



(31) 



lim 9 = 7T , lim ,4 = (32) 

s->0 



while for the boundary S = R/k we choose, 



limi = ~ (33) 
lim 6» = 6» or lim 0' = 9' (R) , (34) 

The field equations ([29]) and ([30]) are the equations presented in (17) and (18) 
of Ref.|5j. The solutions to these equations can be analyzed numerically. From the 
numerical point of view, even if one were interested in the infinite plane, the equations 
have to be solved in a finite disk and eventually analyse the behavior for R — > oo . 



6 



There are basically two methods to solve numerically a boundary value problem 
for a system of ordinary differential equation, generally referred as the "shooting" and 
"relaxation" methods (see Ref. [9]). The first one is easier to implement but has more 
limitations. It consists essentially in guessing the conditions at one boundary to repro- 
duce, after a standard numerical integration (like Euler or Runge-Kutta) , the expected 
behavior at the other boundary. This method is nevertheless of limited usefulness in 
problems with exponentially divergent solutions (as the problem at hand). This is 
so because an infinitesimal error at the initial condition exponentiates to the wrong 
behavior at the other end. Thus, there exists a maximum radius R (which depends on 
machine precision) that can be explored using this approach. The authors of Ref. [5] 
have used a variation of this method to study these equations. 

Indeed, the method implemented in Ref. [5], works as follows. One first takes 

6> (0) (s) = 2 arctan(-) (35) 

s 

to replace at rhs of Eq (|29p . and solves for A(s), choosing the value ^4'(°)(0) in order to 
reproduce the the expected behavior at s = S. With this solution, called A^(s) one 
goes to Eq (|30p to solve for 9(s) giving a solutions that we call O^^s). This iterations 
is carried over on many steps, obtaining solutions, A^(s), 9^(s) in each case one need 
to adjust the value of A'^\0), 8'^{0). The advantage of this method over a standard 
"shooting" approach is that here one needs to guess one initial condition per iteration 
instead of trying to guess two boundary conditions. As mentioned before, due to the 
presence of exponentially divergent solutions at infinity, "very fine tuning" in A'^(0) is 
needed to reproduce the expected —N/ s behavior at "infinity" . In their paper, authors 
quote results for N = 1, S = 30. 

The same boundary value problem can be analyzed more efficiently using a "re- 
laxation" method (see for instance Ref. [9] for details about its implementation). In 
this case, arbitrary disk sizes can be explored since the routine works trying to achieve 
global convergence starting from an initial guess for the solutions in the full interval. 

We have solved the problem using both methods, the one proposed in Ref. [5] and 
the relaxation method mentioned before, for N = 1 and S = 30, imposing as boundary 
condition a(S) = —jj,6(S) = 0. 
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Figure 1: Profile of the gauge field A(s) (left panel) and 9(s) (right panel) for disk size 
S = 30 and N = I. 

In Fig [1] we show the solutions for the gauge field obtained by both methods. After 
30 iterations of the relaxation method, the agreement between the solutions obtained 
by both methods is excellent and are indistinguishable in the graph. They also look 
compatible with the results presented in Ref. [5j. 

As mentioned before, the shooting method becomes inefficient for larger radius 
while the relaxation method we can be use to analyse the behavior of the solutions as a 
function of disk size S. In Fig [2] we show the behavior of the solutions for different disk 
sizes. This figure clearly shows that as the disk sizes increases the solution tends to the 
trivial solutions. Our results clearly contradict the conclusion of Ref. [5] that claim 
the existence of non trivial stable finite energy Chern-Simons CP(1) solitons in R 2 . 
Indeed, our result points in the direction that conclusion of Ref. [3], where instability 
was shown at the perturbative level can be extended to the non perturbative case too. 

The fact that the solutions delocalise as S is increased could be related to the fact 
that the energy decreases with size. In FigJ3]we show the behavior of the energy as a 
function of the disk size, for N = 1. This result can be shown also analytically. Indeed, 
consider the following configuration defined in the interval of length \S 
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s 

Figure 2: (Color online)' Magnetic field profile as a function of scaled radial coordinate s 
for different disk sizes, from top to bottom, S = 30, 60, 120, 240, 480, 960 
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Figure 3: The energy of the solution as a function of disk size S, for N = 1 
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functional (|3ip for the configuration (J36J) in an interval of length XS 

E(XS) = 2tt£ S sds( ^i^fl + d s A xs (s)\ + l{d„Q) 
, (n 2 , 2iVi A5 ( s )\ a £ , ~ 2 



i A s(*) = ^p, 0As(s)=%(|) (36) 

Here A is a real number such that A > 1 and the configurations (I36|) satisfy the boundary 
conditions A\s(s) = — , 9\s( s ) = We can evaluate the energy 

S— Ao s— AS 1 



2 



\ ** + 8 ) C ° S ( 2 } + Axs{s) ) (37) 

We denote the solution corresponding to the interval \S as A\s(s) and 6\s(s) and its 
energy as E(XS). Since the configuration (|36p satisfy the same boundary condition as 
A\s(s) and 0\s(s), we have that 

E(XS) < E(XS) (38) 
Under the transformation s = xX the functional (|37|) becomes 

£(A£) = 2<kJ* xdx{^ f^M + d x A s (x)\ 2 + - A (d x ef 

( N 2 2NA s (x)\ 9 ,0. 9 . ,\ 
+ + x) ( 2 } + ^ (X) ) (39) 

The last expression differ from E(S) only on the factor therefore 

E(XS) < E(S) (40) 

Comparing (|38p and (|40p we have that 

E(XS) < E(S) (41) 

that is, the energy decreases when we enlarge the interval S. 

As we mentioned before, in a finite disk there is some freedom on the boundary 
conditions imposed on the fields. The conditions we have been using below, 6(S) = is 
such the current J<j>(S) = (see Eq. [T7J). The behavior of the energy for more general 
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Figure 4: The energy as a function of the current at the boundary, for S = 30 and N = 1. 

conditions can be easily explored. We show in Fig (H the energy as a function of the 
current at the boundary. 

We would like to mention that the solutions are stabilized if additional potential 
term as in [ID], [11] are included. In particular, we have checked the results of reference 

using the relaxation method and we have found results in complete agreement with 
them. 

In this note then analyzed in more detail the properties of non trivial classical 
solutions in the 2 + 1 CP(l)-Chern-Simons model. Contrary to previous claims, we 
have shown that these solutions are not stable in the infinite plane case, thus, extending 
to the non perturbative domain the conclusions of Ref. [4]. The solutions are perfectly 
well defined in a hnite disk and this could be the case of interest specially if thinking in 
condensed matter applications. The behavior of the solutions under different external 
conditions can be easily obtained following the methods described in this work. 
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